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Abstract 

In this study, we demonstrate closed-loop motion control of self-propelled microjets under the influence of external 
magnetic fields. We control the orientation of the microjets using external magnetic torque, whereas the linear motion 
towards a reference position is accomplished by the thrust and pulling magnetic forces generated by the ejecting oxygen 
bubbles and field gradients, respectively. The magnetic dipole moment of the microjets is characterized using the U-tum 
technique, and its average is calculated to be 1.3 x 10~^° A.m^ at magnetic field and linear velocity of 2 mT and 100 |im/s, 
respectively. The characterized magnetic dipole moment is used in the realization of the magnetic force-current map of the 
microjets. This map in turn is used for the design of a closed-loop control system that does not depend on the exact 
dynamical model of the microjets and the accurate knowledge of the parameters of the magnetic system. The motion 
control characteristics in the transient- and steady-states depend on the concentration of the surrounding fluid (hydrogen 
peroxide solution) and the strength of the applied magnetic field. Our control system allows us to position microjets at an 
average velocity of 1 15 |im/s, and within an average region-of-convergence of 365 |im. 
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Introduction 

Several types of catalytic [1-5] and magnetic [6,7] microrobots 
have been demonstrated to overcome Brownian motion at low- 
Reynolds number regimes [8,9] . Motion of these microrobots is 
based on several propulsion mechanisms. Some of these mecha- 
nisms depend on the efficient transformation of chemical energy 
into kinetic energy using catalytic reaction. Gibbs et al. [10] 
presented a model to explain the driving force for spherical 
colloids, and further showed that the propelling force depends on 
the surface tension of the liquid and the bulk concentration of 
hydrogen peroxide. Solovev et al. [1 1] and Mei et al. [12] have put 
forward a propulsion mechanism based on the catalytic decom- 
position of hydrogen peroxide by microtubular layers of platinum 
and silver, respectively. It has been also shown that self-propelled 
microjets can selectively transport relatively large amounts of 
particles on-a-chip and Murine Cath.a-differentiated cells by 
controlling the magnetic fields in open-loop [13]. Some propulsion 
mechanisms are based on pulling the magnetic microrobots with 
the magnetic field gradients [14-18]. However, these mechanisms 
are limited by the projection distance of the magnetic field 
gradients. Another propulsion mechanism was proposed by Bell et 
al. [19] and Sendoh et al. [20]. This mechanism mimics the 
morphology of the bacterial flagellum and capitalizes on the 
generation of rotating magnetic fields. The rotating fields allow the 
helical microrobot to rotate and act like a corkscrew [8], and 
hence moves without pulling by the magnetic field gradients. 



Biological microrobots [21] have been used by Martel et al. to 
execute non-trivial tasks such as micro-manipulation [22], micro- 
assembly [23], and micro-actuation [24]. The propulsion mech- 
anism of the biological microrobots depends on rotating their 
helical flagella to move forward or backward along the external 
magnetic field lines [25]. This propulsion mechanism allows 
biological microrobots to benefit from the larger projection 
distance of the magnetic field. However, the propulsion force of 
these microrobots is relatively small (10~^^ N) [26], and this 
decreases their range of applications. Wang et al. [27] has studied 
the use of both fuel-driven and fuel-free, as well as ultrasound- 
driven propulsion mechanisms at nano- and micro-scales. 
Although the previously mentioned propulsion mechanisms 
provide solutions to key challenges in microrobotic applications, 
progress towards practical implementation has not yet been 
accomplished. Substantial efforts have been dedicated to develop 
and characterize new propulsion mechanisms, whereas only a few 
attempts have been reported to realize practical applications based 
on closed-loop (feedback) control systems [28]. 

In this work, we demonstrate for the first time the closed-loop 
motion control of self-propelled microjets under the influence of 
controlled magnetic fields generated by a magnetic-based manip- 
ulation system. Figure 1 shows a control exerted over a microjet 
moving towards a reference position under the influence of the 
self-propulsion force and the applied magnetic fields. The ejecting 
oxygen bubbles provide thrust force which allows for the 
navigation of the microjet along the external magnetic field lines 
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(blue lines). We devise a sliding-mode control system [29] owing to 
its robustness in the presence of parameter uncertainties and 
unmodeled disturbance forces such as wall and surface effects, 
bubbles-microjet interactions, and microjet-microjet interactions. 
First, the magnetic dipole moment is characterized based on the 
motion analysis of the microjets using uniform magnetic field 
reversals [30,31]. Second, the characterized magnetic dipole 
moment is used in the realization of a magnetic force-current 
map of the microjet. The magnetic-based control system 
capitalizes on this map. Characteristics of our control system are 
evaluated in the transient- and steady-states using the average 
velocity and average region-of-convergence of the controlled 
microjets, respectively. In addition, the control characteristics of 
the sliding-mode control are compared to the control character- 
istics of a proportional-derivative (PD) control system. The novelty 
of our work lies in demonstrating precise point-to-point closed- 
loop control of self-propelled microjets using weak magnetic fields 
(2 mT). 

Closed-Loop Motion Control of Self-Propelled 
Microjets 

Self-propelled microjets are fabricated by roUed-up nanotech 
[4,32] using layers of platinum, titanium, and iron. These 
microjets are immersed in a solution of hydrogen peroxide 
(H2O2) and surrounded by an array of independent electromag- 
netic coils (Figure 1). The inner platinum layer allows for the 
catalytic break-down of the hydrogen peroxide into oxygen and 
water. A thrust force is generated upon the accumulation and 
release of the oxygen bubbles from one end of the microjet [5]. 
This force allows the microjet to overcome the Brownian diffusion 
and the viscous drag forces. The iron layers allow the microjet to 
orient itself along the external magnetic field lines using magnetic 
torque. Figure 2 demonstrates the motion of a microjet under the 
influence of controlled magnetic fields in open-loop. Throughout 



the experimental work of this study, hydrogen peroxide solution 
with concentration of 5% was used, along with small amounts of 
isopropanol and Triton X. A microjet with a length of 50 |im is 
immersed in the hydrogen peroxide. Uniform magnetic fields (field 
strength of 2 mT at the center of the electromagnetic array) are 
generated using two active electromagnetic coils at a time, and 
used to achieve a square trajectory of the microjet, as shown in 
Figure 2. Magnetic fields are measured using a calibrated three- 
axis Hall magnetometer (Sentron AG, Digital TeslaMeter 3MS1- 
A2D3-2-2T, Switzerland). 

Characterization and Control Law 

In order to position the self-propelled microjet within the 
vicinity of a reference point, a closed-loop motion control system is 
designed (Figure 3). The nominal equation of motion of the 
microjet is given by 



F(P) + f(P,0 + Fci(P) = 0, 



(1) 



where F(P), f(P,05 and Fci(P) are the magnetic force, self- 
propulsion force, and drag force of the microjet at point (P), 
respectively. Our motion control strategy is based on controlling 
the magnetic force (F(P)) towards a reference position. This 
control allows the field lines to be directed towards the reference 
position. A microjet aligns itself along the field lines and moves 
towards the reference position using its self-propulsion force. The 
magnetic force is generated using an external magnetic field 
(B(P)), and is given by [16,26] 



F(P) = (m V)B(P) = (m V)B(P)I = A(m,P)I, 



(2) 



where B(P) is a matrix that maps the input current (I) onto 
magnetic field (B(P)). Further, A(m,P) maps the input current 




50 um 



Figure 1. Closed-loop motion control of a self-propelled microjet under the influence of the controlled magnetic fields. Magnetic 
fields are generated using the magnetic-based manipulation system, shown in the inset at the bottom left corner [14]. The magnetic fields are used 
to orient the microjet towards the reference position (black circle). The microjet moves towards the reference position using the thrust force 
generated by the ejecting oxygen bubbles from its end. The velocity of the microjet in this representative experiment is 100 jim/s for hydrogen 
peroxide solution with concentration of 5%. The inset at the upper right corner shows a Scanning Electron Microscopy image of a microjet fixed to its 
substrate. 

doi:1 0.1 371 /journal.pone.0083053.g001 
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Figure 2. Open-loop motion control of a self-propelled microjet using uniform magnetic fields. The microjet follows a square trajectory 
with an edge length of 500 \irr\, at a linear velocity of 1 00 \xrr\/s. Uniform magnetic fields are generated using 2 active electromagnets at a time. The 
bottom row provides finite element (FE) simulations of the uniform magnetic fields utilized in this open-loop control experiment. // for (/=!,... ,4), 
denotes the current at each of the electromagnets. The magnetic field strength within the center of the workspace is 2 mT. The black arrows point at 
the microjet. 

doi:10.1371/journal.pone.0083053.g002 



onto magnetic force. This force-current map depends on the 
magnetic dipole moment (m) of the microjet and its position (P). 
We estimate the magnetic dipole moment of the microjets to 
realize the magnetic force-current map (2). A microjet aligns itself 
along the magnetic field lines. Reversing the direction of the 
magnetic fields causes the microjet to follow a ^7-turn trajectory, as 
shown in Figure 4. The diameter (D) of the t/-turn trajectory is 
given by [30,31] 



D- 



a7r|P| 
|m||B(P)|' 



(3) 



where a is the rotational drag coefficient of the microjet, and P is 
its linear velocity. In order to determine the magnetic dipole 
moment of the microjet using (3), uniform magnetic fields are 
applied using electromagnets A and G (Figure 1). These fields are 
reversed to initiate t/- turns of the microjets. We repeated this 
experiment 10 times and the average ^7-turn diameter is 



determined from the motion analysis of the microjets. The 
average magnetic dipole moment is calculated to be 
1.31 X 10 ^^A.m^, at magnetic field and linear velocity of 2 mT 
and 100 |im/s, respectively. In (3), a is calculated using [33] 



nf]P 



In 



-0.92 



-0.662 



(4) 



where rj is the dynamic viscosity of the hydrogen peroxide solution. 
Further, / and d are the length (50 |im) and diameter (5 |im) of the 
microjet, respectively. 

The magnetic dipole moment is used in the realization of the 
magnetic force current map (2), this map is a basis of magnetic- 
based closed-loop control systems. First, we devise the following 
sliding-mode control input [29]: 



F(P)= -f(P,0-o();Pref-o(»?r(P), 



(5) 
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Figure 3. Closed-loop control system for precise positioning of self-propelled microjets (inset). P and Pref denote the position of the 

microjet and the reference position, respectively, and I denotes the current vector. 

doi:10.1371/journal.pone.0083053.g003 
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Figure 4. Characterization of the magnetic dipole moment of the self-propelled microjet (blue arrow) using the 6^turn technique. 

The microjet follows a L/-turn trajectory when the magnetic field is reversed. Diameter of the L/-turn trajectory is used to determine the magnetic 
dipole moment. In this representative experiment the diameter is 80 |im, and the magnetic dipole moment is calculated to be 1.48 x 10~^°A.m^ at 
magnetic field and linear velocity of 2 mT and 100 \xrr\/s, respectively. The average magnetic dipole moment is 1.31 x 10~^°A.m^, which is deduced 
from 10 different L/-turn trials using (3) and used in the realization of the magnetic force-current map (2). The red lines represent the magnetic field 
lines. The inset shows a Scanning Electron Microscopy image of a microjet fixed to its substrate. 
doi:10.1371/journal.pone.0083053.g004 

where f(P,0 is an estimate of the propulsion force given in (1). 
Further, Pj-ef is the reference velocity, and r(P) is an additional 
force input and is given by 

r(P)=i(-S(P,0-/?)sgn(s(0). (6) 

In (6), s(0 is the sliding-line and is given by 

s(0 = e + Xe = 0, (7) 

where e and e are the position and velocity tracking errors, 
respectively. Further, and X are positive gains. Finally, ^ (P,0 is 
given by 



/=15s t = 20s t = 22s 

o o 

100 ^im 

^ = 28s (7) / = 33s / = 35s 

25 



Figure 5. Closed-loop control of a self-propelled microjet using the sliding-mode control law (5). The microjet moves towards a 
reference position (small blue circle) under the influence of the self-propulsion force and the magnetic fields. The generated magnetic fields using our 
sliding-mode control system orients the microjet along the magnetic field lines. The microjet moves along these lines using its propulsion force. In 
this representative experiment, the control system positions the microjet at an average velocity of 90 jim/s, and within a region-of-convergence of 
260 jam in diameter. The controller gains are: X=12 and P=13. The large blue circle is assigned by our feature tracking software [14] and the red line 
represents the velocity vector of the microjet. The inset shows a Scanning Electron Microscopy image of a microjet fixed to its substrate. 
doi:1 0.1 371 /journal. pone.0083053.g005 



Motion Control Results 

The sliding-mode control law (5) is implemented by setting (5) 
equal to (2), and solving for the current vector (I) using the 
pseudoinverse of A(m,P). Figure 5 shows a representative closed- 
loop motion control result of the control law (5). The position error 
is determined by tracking the motion of the microjet (large blue 
circle in Figure 5). This position tracking error along with its time 
derivative are used to calculate s(t). Further, the self-propulsion 
force estimate (f(P,t)) is modeled using a periodic force function 
with a period of 0.08 second. This period is determined from the 
motion analysis of the microjets. This time is based on the average 

r = 23s / = 27s ^ 

bp o 
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Figure 6. Representative closed-loop control of a self-propelled microjet under the influence of the controlled magnetic fields, 
using the sliding-mode control law (5). (a) The sliding-mode control system positions the microjet at an average velocity of 97 [irr\/s, and within 
a maximum region-of-convergence (ROC) of 524 \irr\ in diameter. ROCs are represented using the red-dashed circles. The controller gains are: X=12 
and P=13. The red arrows indicate the direction of the microjet, whereas the blue arrows indicate the reference positions, (b) Motion of the microjet 
towards three reference positions (blue vertical lines) along x-axis. (c) Motion of the microjet along y-ax\s. 
doi:10.1371/journal.pone.0083053.g006 



elapsed-time of the ejecting oxygen bubbles. Control law (5) allows 
the magnetic field to orient towards the reference positions (small 
blue circle). The microjet moves towards the reference position at 
an average velocity of 90 |im/s. Due to the self-propulsion force of 
the microjet, our control system cannot achieve zero position 
tracking error. It rather positions the microjet within the vicinity of 
the reference position, which we refer to as a region-of- 
convergence as shown in Figure 6. 

Figure 6 shows the controlled motion of a microjet towards 
three reference positions. Within the vicinity of the reference, the 
control system reverses the direction of the magnetic field based on 
the position and velocity tracking errors (7). This magnetic field 
reversal allows the microjet to stay within the vicinity of the 
reference position. Diameter of this region-of-convergence de- 
pends on the linear velocity of the microjet, the magnetic torque 
exerted on its magnetic layers, and the bandwidth of the control 
system (our microscopic vision system has a maximum frame-per- 
second of 15). The sliding-mode control gains {X and j5) must be 
positive. The gain X controls the slope of the sliding-line (7), 
whereas the gain controls the convergence time of the errors (e 
and e) to the sliding-line. In this representative experiment 



(Figure 6), we observe that setting the controller gains to X=\2 
and j^=1.3, results in an average velocity of 97 |im/s and 
maximum region-of-convergence of 524 |im. This motion control 
trail is repeated 10 times, and the average velocity and average 
region-of-convergence are calculated to be 115 + 26 |im/s and 
356+ 110 |im, respectively. 

In order to evaluate the characteristics of the sliding-mode 
control system, we compare its results to a conventional PD 
control system. In this case, the magnetic force is given by 

F(P) = Kpe + Kde. (9) 

where Kp and Kd are the controller positive-definite gain matrices. 
Control law (9) is implemented by setting the PD control force 
equal to the magnetic force-current map (2), and calculating the 
current vector using the pseudoinverse of A(m,P) that is based on 
the characterized magnetic dipole moment and position of the 
microjet. We observe that the PD control system achieves an 
average velocity and region-of-convergence of 1 19 + 30 |im/s and 
417+115 |im, respectively. The averages are calculated from 1 0 
closed-loop motion control trials. Figure 7 shows a representative 




^ fnir"! .V [mm] y [mm] 

(a) Motion of a controlled microjet (b) Motion along .v-axis and time (c) Motion along i -axis and time 



Figure 7. Representative closed-loop control of a self-propelled microjet under the influence of the controlled magnetic fields, 
using the proportional-derivative (PD) control law (9). (a) The PD control system positions tlie microjet at an average velocity of 127 jim/s, 
and within a maximum region-of-convergence (ROC) of 662 \irr\ in diameter. ROCs are represented using the red-dased circles. The controller gains 
are: kpi =/cp2 = 15 and k^i = k^2 = 5, where {k^i) and {k^t) are the entries of the gain matrices (Kp) and (Kd) for (/= 1,2), respectively. The red arrows 
indicate the direction of the microjet, whereas the blue arrows indicate the reference positions, (b) Motion of the microjet towards three reference 
positions (blue vertical lines) along x-axis. (c) Motion of the microjet along y-axis. 
doi:1 0.1 371 /journal.pone.0083053.g007 
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motion control result of the microjet using the PD control system. 
These results show that the non-linear sliding-mode control system 
achieves 14% higher positioning accuracy, as opposed to the linear 
PD control system, in the steady-state. 

Conclusions 

Closed-loop control of self-propelled microjets is implemented 
using a magnetic-based manipulation system and a sliding-mode 
controller. Magnetic fields are generated and oriented based on 
the position tracking error of the microjet with respect to the 
reference position. This control allows the microjet to orient along 
the field lines using the magnetic torque exerted on its iron layers. 
The magnetic torque is calculated to be 2.6 x 10~^^ N.m, based 
on the characterized magnetic dipole moment of the microjet and 
the maximum magnetic fields (2 mT) used throughout the 
experimental work. This magnetic torque allows the microjet to 
overcome a drag torque of 3.3 x 10 N.m, based on an angular 



velocity of 25 rad/sec. Our control system uses the magnetic 
torque only to keep the microjet within the vicinity of the reference 
position due to the self-propulsion force of the microjet that cannot 
be controlled during a motion control task. The unmodeled 
dynamics due to the wall and surface effects and microjet-bubbles 
interactions are mitigated using a sliding-mode control system. 
This control system achieves 14% higher positioning accuracy 
than a magnetic-based PD control system. The accuracy provided 
by the proposed closed-control system could allow microjets to be 
used in micro-manipulation, micro-assembly, micro-actuation, 
and applications that are not yet conceived. 
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